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A model for the study of the effective quasistatic conductivity and permittivity of dispersed systems with particle- 
host interphase, within which many-particle polarization and correlation contributions are effectively incorpo- 
rated, is presented. The structure of the system's components, including the interphase, is taken into account 
through modelling their low-frequency complex permittivity profiles. The model describes, among other things, 
a percolation-type behavior of the effective conductivity, accompanied by a considerable increase in the real 
part of the effective complex permittivity. The percolation threshold location is determined mainly by the thick- 
ness of the interphase. The "double" percolation effect is predicted. The results are contrasted with experiment. 
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1. Introduction 

The studies of disperse systems, nanofluids, and systems of nanoparticles are nontrivial, but impor- 
tant for many fields in science and industry. Much attention has been recently focused on the question 
how the properties of the particle-host interphase affect the dielectric characteristics of the entire sys- 
temjll|; 40. This is an essentially many-particle problem, playing a crucial role, say, in percolation stud- 
ies dfoi. The latter involve (a) various numerical (mainly Monte-Carlo) methods 00], (b) graph jioll and 
renormalization group 1 7] theories, (c) different improvements of the Maxwell-Garnett (ll|l and Brugge- 
man 1 12] approaches. 

It should be noted that the practical realization of methods (a) is labor-consuming and imposes high 
computer specifications. Methods (b) provide certain scaling relations, but ignore the physical nature 
of dielectric response of an inhomogeneous system. Methods (c) seem more adequate physically, but 
are actually one-particle approaches. They can demonstrate a percolation-type behavior of conductiv- 
ity but ignore the fractal nature of percolation cluster flill . Numerous attempts at improving these 
approaches are widespread in modern technical literature (see Mill ), but are, for most part, poorly justi- 
fied physically. 

In the present work, the role of the interphase in the formation of a percolation-type behavior of 
conductivity is analyzed within an electrodynamic model based upon the method of compact groups of 
inhomogeneities I16H18I1 . The "compact group" means any macroscopic region, in which all distances 
between the inhomogeneities are much shorter than the wavelength of probing radiation in the medium. 
The method allows one to effectively sum up and average the electric field and displacement vectors 
in finely-dispersed systems of particles with complex internal structure and arbitrary permittivities. In 
doing so, unnecessary detailed elaboration of the processes of interparticle polarization and correlations 
in the system is avoided. 

The main idea of the method is that in the long-wave limit, a macroscopically homogeneous and 
isotropic particulate system is electrodynamically equivalent to an aggregate of compact groups of par- 
ticles. The compact groups include macroscopic numbers of particles and can be viewed as one-point 
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inhomogeneities, with negligibly small particle number fluctuations inside. They are in fact self-similar. 
Note that this property is characteristic of a percolation cluster. 

The model of the disperse system under consid- 
eration is depicted in figure [T] Each particle con- 
sists of a spherical hard core (black), with diame- 
ter d and permittivity E\. The core is surrounded 
by a penetrable concentric shell (gray), with outer 
diameter d{\ + 8) (<5 being the relative thickness of 
the shell) and permittivity £2- The particles are em- 
bedded into a host medium (white), with permittiv- 
ity £q. All the permittivities are in general complex. 
The volume concentration of the hard-cores is c. The 
shells can freely overlap; together, they form the in- 
terphase between the hard cores and the host. The 
volume concentration <p- mt of the interphase for par- 
ticular values c and S can be calculated by statisti- 
cal methods. We take it from Monte Carlo simula- 
tions in, Hob. 

We expect that, in the main features, our results 
should remain valid for macroscopically homoge- 
neous and isotropic systems with non-spherical 
core-shell particles: the sphericity comes into play 
at the points where the explicit form of the func- 
tion (pi nt = </>i nt (c,5) is needed. By contrast, spherical 
symmetry of the particles is crucial for most approaches to the problem which often represent different 
versions of one- and two-particle approximations and use it either implicitly or through the require- 
ment that the shells should be hard. On the other hand, the systems of anisotropic particles are usually 
treated as if the latter were solitary. Some examples and further literature on these topics can be found 
in II El El El. 




Figure 1. The model under consideration. 



2. General equation for the low-frequency complex permittivity 

The structure of the system's components, including the interphase, is taken into account through 
modelling their complex permittivity profiles and the choice of the way of homogenization. 

We are interested in the low-frequency value e = e' — \Analu> of the effective complex permittivity 
and determine it as the proportionality coefficient in the relation 

<D(r)> = <e(r)E(r)> = £-<E(r)>, (2.1) 

where D(r), E(r), and e(r) are the local values of, respectively, the electric displacement, electric field, and 
permittivity in the system, and the angular brackets stand for the statistical averaging. In the present 
work, we suggest that this averaging can be reduced to the following two-step procedure. 

1. Finding the equilibrium distribution of the particles with penetrable outer shells. This distribu- 
tion is characterized by the hard-cores' volume concentration c and the particles' effective vol- 
ume concentration (p - (p{c,5) (the latter is the sum of c and the volume concentration of the re- 
gions occupied by the shells). Electrodynamically, such a system can be viewed as an aggregate of 
non-overlapping (white, black, and gray in figure [1} regions with permittivities £; = e^. - \AnOi Id), 
i = 0, 1,2, and volume concentrations l-cp, c, and (p- c, respectively. 

2. Electrodynamic homogenization, by direct integration with respect to its volume of this ag- 
gregate within the compact group approach. This suggests that: (a) with respect to a probing field 
of long wavelength A, the aggregate itself is viewed as a set of regions (termed as compact groups) 
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with typical linear sizes much smaller than A, but yet remaining macroscopic; (b) the macroscopic 
dielectric properties of the aggregate are equivalent to those of a macroscopically homogeneous 
and isotropic system prepared by embedding the aggregate's non-overlapping regions into some 
fictitious medium, of permittivity e. The permittivity distribution in this fictitious system can be 
written as 

e(r) = e + <5c(r), (2.2) 

where, in our case, 

N N' 

Se(r) = (e„ - e) n (r, Q ) + E ( £ i " £ ) n i ( r > n a) + L - <0 U 2 (r, n' b ) (2.3) 

a=l b=\ 

is the local permittivity deviation from e due to the presence of a compact group of the aggregate's 
non-overlapping regions at point r. Here, Flo (r,Qrj), 111 (r,Q a ), and 1I2 (r.n^) are the characteristic 
functions of the regions Qrj> and D.' b occupied by, respectively, the real host, the ath particle's 
hard core, and the Mh connected cluster of the overlapping outer shells: 

,' 1 , r £ O, 

n,-(r,Q) = ^ (2.4) 
0, rCQ. 



With the model permittivity distribution <2.2> and (273) proposed, the averaged displacement and field 
in formula d2.lt are calculated by the general formulae 



s=l 



(E) = \l+^\-^\ <[5c(r)] s >[E , (2.5) 



(D> = |e + e£[-^] ([de{r)] s ) + p [-^j <[5c(r)] s+1 > J E . (2.6) 
In view of the property <2.4t , the averages 

<[5e(r)] n > = i f [5e(r)]"dr 

v 

are easy to find: 

<[<5e(r)]"> = (1 - </>)(£„ - e) n + c{e 1 -e) n + {<p-c) (e 2 -e) n . (2.7) 
Substitution of i2.7) into series <2.S) . <2.6) and summation of those give 



£n-e ei— e £?-£ E — e 

[l-d))— + c— + (ch-c)— = . (2.8) 

2e + £ 2e + £i 2e + £ 2 2e + £ 

The permittivity e is the only unknown parameter left. In what follows, we suggest that the dielec- 
tric properties of the fictitious and effective systems are identical, that is, e = £. This approximation is 
equivalent to the Bruggeman-type homogenization, when the host and dispersed particles are treated 
symmetrically. This seems reasonable for the system under study because the shape of the white and 
gray regions is extremely complicated. The degree of this complexity increases with c even more. 

The other extreme case, e = £q, corresponds to the Maxwell-Garnett type of homogenization. 



3. Main features of the model 

In the quasistatic limit and at e = £, the general equation (2.8) reduces to two real equations, a cubic 
one for a and, once a is found, a linear for e'\ 

(Jq — ct <j\—(j ay — <J 
(1-0)— + c— + {(p-c)— = 0. (3.1) 

2(7+ (Tq 2(7 + <7 1 2CT + CT2 
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a-(p)— t + c- t + (0-c)— 7 = 0. (3.2) 

(2cr + cro) 2 (2<r + (Ti) 2 {2a + o 2 ) 2 



A general analytical analysis of equation <3.1> can be carried out with the use of Cardano's formu- 
las and is cumbersome jgjfl . Nonetheless, the major features of the model can be grasped with simple 
physical reasoning. For the sake of convenience, we change in equations (373} and (3.2) to dimensionless 
variables x = alo\, y - e 1 ' le' and x\ = Oilo\, yi - e'^Eq U = 0, 1,2). Typically, xq « 1 and y > 1. 



3.1. Percolation-type behavior 

In the limit of a non-conducting matrix, xq — ► 0, equation d3.1> has three solutions, x = and 




(3.3) 



If X2 > 0, a physically meaningful nontrivial solution (that with the positive sign in front of the square 
root, see figure[2) appears only under the condition 

0(c c ,5) = i (3.4) 

and is independent of X2. In the case X2 = 0, the well-known value of c c = 1/3 for the two-component 
Bruggeman model occurs fiUl . 
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Figure 2. Percolation (dashed line, 5 = 0) and "double" percolation (solid line, S = 0.05); xq = 1 • 10 
x 2 = 5-10" 5 . 



-10 



The relation (3.41 determines the threshold concentration c c for the effective conductivity. This value 
is affected only by the fact of the existence of the interphase (8 + 0), rather than its dielectric and conduc- 
tive properties. 

Our estimate of c c as a function of the relative thickness 8 of the interphase layer is shown in figure[3] 
The calculations were performed using the Monte Carlo results 12011 (A; = (1 + 8)~ 3 , q> = elk): 

(P = 1 - (1 - lap) exp [-(1 - k)cp] ex P |- 2(1 fcy fcy)3 [( 8 ~ 9fcU3 + fc ) 

- (4+ 9fc 1/3 - 18fc 2/3 + 5k)k<p + 2(1 - k)k 2 <p 2 ] \. (3.5) 
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Figure 3. Threshold concentration as a function of the relative thickness of the interphase (solid line). 



Dashed line: c c = i (1 + 8) 3 



The analysis revealed that for realistic c < 0.5, the relation c c = | (1 + 6) 3 can be used. 

In the immediate vicinity of c c (c — ► c c + 0) and for nonzero 8, formula d3.3> takes the form 



4 



1 + 



i + c(l-x 2 ) 
■c(l-x 2 ) 



1 

y 3 



(3.6) 



Correspondingly, the effective conductivity o oc (c- c c ) f where the critical exponent t a 1.0. The effective 
permittivity e', as follows from equation <3.2h increases anomalously as xq — » 0. The latter fact is in accord 
with predictions Esll . 



3.2. Effective critical exponent of conductivity 

In practice, both the threshold concentration c c and the critical exponent t are determined by inter- 
polating the conductivity experimental data a - a[c), obtained for some finite interval c e [c\,C2] near c c 
(ci > c c ), with the scaling law a = A(c - c c ) 1 , A and t being independent of c. Then, 



, c(c 2 ) /, c z -c c 



(3.7) 



and this value is expected to be a c-independent constant. 

By contrast, the asymptotic formula <3.6> reveals that even small variations of c near c c cause the 
expression in the brackets, which is proportional to A, to change considerably. This means that a formal 
application of the above procedure and formula H3.7t to a system with conductivity 13.3) will result in an 
effective exponent f e ff sensitive to the parameters c\ and c 2 (figure[4}. In particular, for a given 3^0, f e g 
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Figure 4. Effective critical exponent of conductivity as a function of c 2 at fixed c\, 8 - 0.1 (c c = 0.251), and 
x 2 - 5 • 10~ 5 , calculated with formulas (33} and (TTJ. From bottom to top, c\ = 0.26, 0.27, and 0.28. 
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Figure 5. Effect of the matrix's conductivity on the 
effective permittivity. From top to bottom, xq = 1 • 

a . r . . a — 



10 



1-10 5 , and 1-10 4 . The other parameters: 



yi = 1.5, y 2 = 1, X 2 = 0.05, 8 = 0.05 




Figure 6. Effect of the interphase thickness on the 
effective permittivity. From right to left, 5 = 0, 0.05, 
and 0.10. The other parameters: y\ = 1.5, y 2 = 1, 
x = 1 ■ 10" 5 , x 2 = 0.05. 



increases as the interval [c\,c 2 \ (where c 2 < 1/3) is: (a) shifted to higher values of c (while its width is 
fixed); (b) widened at fixed c\ . Also, the threshold concentration found within this procedure is expected 
to exceed c c . 

Different 3D percolation models d] and renormalization group calculations ^ give for t esti- 
mates of a 1.3 +1.7 and w 1.9, 2.14, respectively. Experimental values of t are usually 1.5 4- 2.0 and some- 
times even twice as much l22Tl . As evident from figure|4] our theory is capable of reproducing a variety of 
these values. 



3.3. Effect of the matrix's conductivity 

For real substances, xq + 0, though xn can be extremely small. As a result, the percolation-type depen- 
dence of x with c changes to a smooth one, with its slope considerably increasing near c c . Simultaneously, 
the maximum value of y becomes bounded above and decreases as xq increases (figure[5). The location 
of the maximum shifts to lower concentrations as S increases (figure[6). Calculations show that it remains 
practically independent of x 2 and actually equal to c c . 

Below the percolation threshold, the effective conductivity is usually approximated by the scaling law 
a = B{c c - c)~ s . Given experimental data for some interval c e [c\,c 2 ], c 2 < c c , the effective values s e ff of 
the critical exponent s are found as 



, cte) /, c c -c 2 
s = "log—— -/log 

C(Ci) ' c c - Cl 



(3.8) 




Figure 7. Effective critical exponent of conductivity below c c as a function of xq for 8 = 0.1 (c c = 0.251) 
and x 2 = 5 • 10~ 5 , calculated with formulas and (3.1) and 13.8) at c\ = 0.24 and c 2 = 0.25. 
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Our estimates of Soff with formulas 13.1) and d3.8t are shown in figure [7] They correlate well with 
typical theoretical j2^|231 and experimental [ 22] values of 0.75 and 0.7 1.0, respectively. 

3.4. Behavior of the permittivity 

According to equation H3.2t . the effective permittivity is given by 

(2x+x ) 2 (2x+Xo) 2 
(1 - d>) e' n + c tt£\ +{(b-c) -E? 

° {2X+l) \ l * X+X f \ (3.9) 

(2x+Xo) (2x+x ) 

(1 - </>) Xo + C + {(p - C) — - -rX Z 

{2x+lY {2x+X2r 

For a badly-conducting matrix (xo ->• 0) and under the condition x « 1, three particular situations are 
of interest to point out to. 

1. The system is below the percolation threshold and the conditions x « ^/xq, x « ^/xqXz, x « xt 
(that is, a « ^oqO\, o « ■ v /Tjo7j2, and a « o 2) are met. Then, the dominant contributions to both 
the numerator and the denominator are made by their first terms, and we expect that e' cx x cx 

2. The system is above the threshold and x » y/xo, x » yfxz, x » X2 (c » v^oci, a » v^i ^. and 
a » <7 2). Now, of significance become the first and the third terms in the numerator, the latter being 
almost independent of x, and the second term in the denominator. Correspondingly, the e' versus 
c dependence is expected to be close to e' cx x -1 cx (c- c c )~ feff , with the proportionality constant 
slightly dependent on c. 

The exponents s e ff and f e ff in the two preceding scaling-like laws are independent of the compo- 
nents' permittivities e' v 

3. The system is close to the percolation threshold, x » ^/xq, and x » X2 (<7 » v /?ooT, a » 02). Then, 
the numerator is almost x-independent, whereas the denominator is mainly contributed to by the 
second and the third terms. The dependence of e' on x takes the form e' cx ax I (l + fox 2 ), where the 
coefficients a and b are easy to recover. 



3.5. Applicability to real systems 

Figure[8]shows the results of processing with formula (379) the experimental data 1 28] for the effective 
permittivity of the composites prepared by embedding spherical Ag particles (a mean radius « 100 A ) into 
a KC1 matrix. The particles were made by evaporating Ag in the presence of argon and oxygen gases so 
as to form a thin (according to the authors, of approximately 10 A , 5 = 0.10) oxide coating on them. This 
coating prevented the particles from cold-welding together, but was thin enough to allow metal-to-metal 
contact under high pressure. The composites were prepared by mixing Ag particles and KC1 powder and 
then compressing the mixture into a solid pellet under high pressure. 

It is seen from figure [8] that formula d3.9t not only reproduces data l28ll over the entire range of Ag 
concentrations investigated, but also gives an estimate of 8 = 0. 14 -h 0.19, sufficiently close to the expected 
one. 

The conductivity (resistivity p) data for a few Ag-KCl composite samples, prepared in the above way, 
are given in Eill . They pertain only to a very narrow vicinity of the percolation threshold, where p drops 
7 orders of magnitude with a 1% Ag volume concentration increment; the parameters of the KC1 matrix 
are not specified. As figure [9] reveals, formula d3.lt can reproduce data 12911 sufficiently well. Better fits 
can be produced by introducing c-dependences for some of the parameters of the model. These facts may 
indicate that, in addition to experimental errors, various other factors and phenomena (inaccuracy of the 
function cf> = <p{c,8), particles' size distribution, silver dissolution and local dielectric breakdown in the 
KC1 matrix, polarization effects, etc.) come into play as c c is approached. The analysis of them goes far 
beyond the scope of this paper. 
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Figure 8. Effective permittivity data [28] for two series of Ag-KCl composite samples (circles and triangles) 
below the percolation threshold and their fits with formula 13.9) at e' q = 5.0, S = 0.186 (solid line) and 
e' Q = 7.0, S = 0.145 (dashed line). The dotted lines are the scaling-type fits (with c c = 0.20, s e ff = 0.72 and 
c c = 0.22, s e ff = 0.74, respectively) proposed in ji^l to the data for c> 0.11. 



3.6. "Double" percolation 



For intermediate values of x 2 (xo « x 2 « x\)> a "double" percolation can be noticeable, that is, a 
new increase in x after some levelling off (figure [2); it is accompanied by the appearance of a new peak 
in the concentration dependence of the permittivity (figure [10). The physical cause of this phenomenon 
is clear — in a concentrated system, the hard cores of particles with penetrable shells begin to contact 
intensively to form percolation cluster and add to the effective conductivity. 

Evidently, the threshold concentration c' c for "double" percolation is close to a value of 1/3. In the 
region |c - 1 1 « X2 « 1, the c-dependence of the effective conductivity H3.31 is represented by the square- 
root law 

1/2 i 



\ (3x 2 ) 



1/2 



+ O(x 2 ) = -[3x 2 0'(c c ,<5)] 



1/2 



Cc- 



c c ) 1/2 + 0(x 2 ), 



(3.10) 



<p' being the derivative of (p with respect to c. For concentrations satisfying the condition c - I » x 2 , it 



logb/ljQ-m)] 




0.21 



0.215 



0.225 



Figure 9. Effective resistivity data 1 29] for Ag-KCl composite samples (squares) and the fit to them with 
formula (33} at o\ = 6.3 • 10 7 S/m, x = 5 ■ 10" 16 , S = 0.162 (c c = 0.214), and x 2 = 4 • 10" 6 . 
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Figure 10. Effective permittivity as "double" percolation occurs; xo = 1-10 8 ,X2 = 5-10 4 ,yi =1.5,y2 = l 
8 = 0.05. 



becomes linear, 




(3.11) 



with a considerably greater amplitude as compared to those in formulas <3.6> and <3.10> . 

As for now, we are unaware of experimental observations of the effect described. Usually, "double" 
percolation is associated with a non-monotonous behavior of the conductivity in composites made by 
embedding the conducting particles into a two-component matrix (see, for example, I30l,l3lll ). 

It can also be shown that a decrease in x can occur at c ~ c' c in systems with xz » 1, Similar effects 
were observed in two-phase composite solid electrolytes with a highly conducting interphase layer [22]. 



The model proposed is interesting in the sense that it is based on rather clear assumptions, incor- 
porates many-particle effects in a consistent way, and can be further refined so as to apply to a more 
complicated systems, including fluctuation phenomena, etc. For practice, it can serve as a rather flexible 
theoretical basis for analysis of dielectric and conductive properties of such complex systems as various 
dispersions, colloids, and so-called nanofluids, or for the development of new composite materials. 



We are grateful to an anonymous Referee for stimulating remarks. 
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npoBiflHicTb Ta flie/iei<TpnHHa npoHUKHicTb flucnepcHnx 

CI/lCTeM i3 BmbHO-npOHMKHMM Mi>l«()a3HMM ujapoM Mi>K 

HacTi/iHKaMi/i ii cepeflOBMineM 

M.fl. CyniKO, A.K. CeMeHOB 

OflecbKww Hai^ioHa^bHwPi yHiBepcwTeT iivieHi 1. 1. Me^HUKOBa, syn. flBopuHCbKa, 2, 65026 Ofleca, YKpaiHa 
3anponoHOBaHO MO,qe/ib fl/ifl BWB4eHHfl e^eKTWBHWx KBa3iCTaTHHHHx npoBiflHOCTi Ta flie/ieKTpw | -iHOi npoHWK- 

HOCTi flUCnepCHHX CWCTeM i3 Bi/lbHO-npOHUKHUM Mi>KC(>a3HHM LUapOM Mi>K HaCTMHKaMU Ta cepeflOBWiueM, flKa 

eefieKTMBHO BpaxoBye 6araT04acTWHKOBi no/iflpn3aLjii/iHi Ta Kope/iflL^MHi ecf>ei<™. CTpyKTypa KOMnoHeHTiB cuc- 

TeMn, BK/lKJ4aK)MW MiX(()a3HWii LUap, BpaXOByETbCfl LU/lflXOM MOfle/llOBaHHfl npO<j>i/liBiX HH3bKOHaCTOTHOI KOMn- 

/ieKCHoTflie/ieKTpn | -iHOi npoHUKHOCTi. Mofle^b, 30i<peMa, onncye nepKoyiai^iwHO-nofli6Hy noBefliHKy ecfieKTHBHOi 

npOBiflHOCTi, LUO CynpOBOflXyETbCfl CVTTEBUM 3pOCTaHHflM flMCHOi 4aCTHHH ecfieKTHBHOi KOMn/ieKCHOI npOHWK- 

HOCTi cwcTeMH. no/ioxeHHfl nopory nepKO/iflL4iT BH3Ha4aeTbCfl to/iobhum hwhom tobluhhokj Mix(|>a3Horo Luapy. 
nepefl6a4eHO ecfieKT "noflBMHoi" nepKO/iflL(iT. Pe3y/ibTa™ nopiBHHHO 3 eKcnepnivieHTOM. 

K/iK>HOBi oiOBa: CTpyKiypoBana ^acmHKa, fli/icnepcHa cticreMa, flieneKTpmHa npoHUKHicrb, npoBiflnicTb, 
nepKO/imj,w 
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